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/******** * * """\ 

* 

* Copyright (c) Schlumberger Technologies 1999. All rights reserved. * 

* No part of this program may be photocopied, reproduced , or * 

* translated to another programming language without the prior * 

* written consent of Yield Enhancement Systems. * 

* 

\********** — * «*" "* * * *********/ 

/******************************** ************************* 

Returns : 

0 for success 

-1 for functional failure 

-2 if algorithm fails to find acceptable threshold ^ 

#include <stdio.h> 
ftinclude <stdlib.h> 
# include <string.h> 
ft include <math.h> 
ft include <mp.h> 
#include < image. h> 
ftinclude <SLBapi.h> 
#include <pgm.h> 
# include <morphl.h> 
# include <scatter2d8 .h> 

ftdefine SQRT2 1.414213562373 /*sqrt(2)*/ 

#define MIN DIST 1 /*min image difference used as the 

~ beginning of ID profiles*/ 

ftdefine MVA_SIZE 21 /*size of moving average window*/ 

ftdefine MVA_RATIO 0.9 /*allowable shrinkage in Id profile 

wrt its mva*/ 

#define MIN_DATA_SIZE 80 /*min ft of data points to use curve 

~ fitting*/ 
ftdefine MX_SIZE (3 << GRAY_SCALE_BITS) /*size of the data matrix for curve 

fitting*/ 

ftdefine HIST_MORPH_PASSES 1 /*no. of passes for open/close 

operations*/ 

ftdefine SIGMA_MAX 5 /*enable removal of outliers if sigma 

of fit is bigger than this*/ 
/*ftdefine MIN THR 10 /*min thr below which no difference 

~ is called a defect*/ 

ftdefine EPSILON 0.1 /*min value for det(AtA)*/ 

ftdefine ALG_FITMASK 1 /*detection alg using cleaned 2D hist 

as a mask*/ 
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static void hist2D8 (unsigned char *a, unsigned long tcols_a / 
unsigned char *b, unsigned long tcolsjb, 
unsigned long ncols, unsigned long nrows, 
unsigned char * table) ; 



static OneD_dist_prof ile distmax, distmin; 

/™ * * * 

int Find2DHistoOutliers(OIP_byte_image *pSrdmagel, 
OIP_byte_image *pSrdmage2, 
OIP_byte_image *pDstIraageResult ( 
PIXEL_diff jpara PixelDif f Para, 
unsigned char *hist, 
unsigned char *hist2DBitMap, 
unsigned long FovID, 
long short test) 

{ 

int rc; 

unsigned char *almagel, *almage2, *aResult; /* ROI pointer */ 

const long Passes « HIST_MORPH_PASSES; /*l=3x3 open operation, 2=5x5, etc*/ 

char name [128] ; 



• 



/* check if the ROIs have the same size */ 

if ( (pSrcImagel->Info.RegWidth != pSrdmage2->Info.RegWidth) 

(pSrdmagel->Info.RegWidth != pDstImageResult->Info.RegWidth) 
(pSrdmagel->Info.RegHeight != pSrcImage2->Info.RegHeight) 
(pSrcImagel->Info.RegHeight != pDstImageResult->Inf o.RegHeight) ) 
return -1; 



/* determine ROI pointers */ 
almagel - pSrcImagel->Buf Ptr 

+ pSrcImagel->Inf o. RegXO 

+ pSrdmagel->Info.RegYO * pSrdmagel->Info. Width; 

almage2 = pSrcImage2->Buf Ptr 

+ pSrcImage2-> Info. RegXO 

+ pSrdmage2->Info.RegY0 * pSrdmage2->Info. Width; 

aResult = pDstImageResult->Buf Ptr 

+ pDstImageResult->Info. RegXO 

+ pDstlmageResult -> Info. Reg YO * pDst ImageResult -> Info. Width ; 

memset(hist, 0, GRAY_LEVELS * GRAY_LEVELS * sizeof (unsigned char) ) ,- 
hist2D8 (almagel, pSrcImagel->Info. Width, 

alraage2, pSrdmage2-> Info. Width, 

pDst ImageResult- >Inf o . RegWidth, 

pDst ImageResult- >Inf o . RegHeight , 

hist) ; 

if (shorttest) 

{ . 

sprintf (name, " . . /data/2Dhisto_%ld.pgm" , FovID) ; 
printf ("Writing 2Dhisto to %s\n" , name) ; f flush (stdout) ; 
rc = write_image (name, hist, 256, 256); 
if (rc) ERROR_RETURN(rc, "write 2Dhisto failed" ) ; 
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return 0; 



void hist2D8 (unsigned char *a, unsigned long tcols_a, 
unsigned char *b, unsigned long tcolsjb, 
unsigned long ncols, unsigned long nrows, 
unsigned char *table) 

/************************************************* 

a,b -- TCL of ROI (images) 

tcols_a -- image width 

ncols, nrows -- ROI size 

table -- scatter table (a 256x256 image) 

****************** i 
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unsigned long i, j; 

unsigned long skip_a, skip_b; 

unsigned long index; 

skip_a = tcols_a - ncols; 
skip_b = tcols_b - ncols; 

for (i = 0; i < nrows ; ++i) 

{ 

for (j = 0; j < ncols; ++ j ) 



index=( (*a++)«GRAY_SCALE_BITS) + (*b++) ; 
table [index) = 1; 



a +- skip_a; 
b += skipjb; 



